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ABSTRACT 

A theory of pulsar radio emission generation, in which the observed waves are produced 
directly by maser-type plasma instabilities operating at the anomalous cyclotron-Cherenkov 
resonance u — k\\v^\ + lub/ Ires — and the Cherenkov-drift resonance oj — k\\v\\ — k±Ud = 0, 
is capable of explaining the main observational characteristics of pulsar radio emission. The 
instabilities are due to the interaction of the fast particles from the primary beam and the tail 
of the distribution with the normal modes of a strongly magnetized one-dimensional electron- 
positron plasma. The waves emitted at these resonances are vacuum-like, electromagnetic waves 
that may leave the magnetosphere directly. In this model, the cyclotron-Cherenkov instability 
is responsible for core emission pattern and the Cherenkov-drift instability produces conal 
emission. The conditions for the development of the cyclotron-Cherenkov instability are satisfied 
for both typical and millisecond pulsars provided that the streaming energy of the bulk plasma 
is not very high 7^ w 10. In a typical pulsar the cyclotron-Cherenkov and Cherenkov-drift 
resonances occur in the outer parts of magnetosphere at r^es ~ lO^cm. This theory can account 
for various aspects of pulsar phenomenology including the morphology of the pulses, their 
polarization properties and their spectral behavior. We propose several observational tests for 
the theory. The most prominent prediction are the high altitudes of the emission region and the 
linear polarization of conal emission in the plane orthogonal to the local osculating plane of the 
magnetic field. 



Subject headings: stars:pulsars-plasmas-waves-radiative transfer 



1. Introduction 

More than twenty five years have passed since the discovery of pulsars and there is still no consensus on 
the basic emission mechanism. At the present time, there are about a dozen competing theories which differ 



both in the physical effects responsible for the radiation and in the locations where they operate ( [Melrose 



1995| ) . Probably the only point of agreement between all these theories is the association of pulsars with 



magnetized, rotating neutron stars. By contrast, there is so much observational data available that none of 
the existing theories can explain all the main observational facts. 

There are several reasons that have precluded understanding of pulsar radio emission. First, there 
is the unusual physical conditions found in pulsar magnetospheres (relativistic electron-positron plasma, 
superstrong magnetic fields, ultrarelativistic beam). Secondly, only a small fraction of the energy lost by 
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a neutron star is re-emitted in the radio (even in the high energy range, where a considerable portion of 



energy is emitted, there is still no consensus on the origin of this emission, e.g., Daugherty & Harding 



1996, Romani 1996). The third reason is that after thirty years of research, we still do not know the 



general structure of a pulsar magnetosphere ( Mestel 1995| ). We understand only particular features, like the 
existence of open and closed field lines and where electric field can be parallel to the magnetic field. 

To date, the most widely discussed theory attributes the emission to coherent curvature emission by 
bunches of particles. Though this theory can explain a broad range of observed pulsar properties by the 
careful arrangement of the magnetic field geometry and of the form and size of bunches, thirty years of 



theoretical efforts have failed to explain the origin of these bunches (Melrose 1995). This theory can also be 
ruled out on the observational grounds (Lesch et al. 1998). In addition to the work of Lesch et al. 1998| we 



note, that this theory also fails to explain the observed correlations of the conal peaks Kazbegi et al. 1991a 
and a large size of the emitting region Gwinn et al.l997 (see Section ||). 



We propose that the pulsar radiation is generated by plasma instabilities developing in the outflowing 
plasma on the open field lines of the pulsar magnetosphere. Plasma can be considered as an active 
medium that can amplify its normal modes. In the case of the two instabilities discussed below, the wave 
amplification is due to the resonant wave-particle interaction, i.e., in the rest frame of the particle the 
frequency of the resonant wave is zero or a multiple of the gyrational frequency. The plasma instabilities 
that we argue operate in the pulsar magnetosphere may be described by the (somewhat contradictory) 
term "incoherent broadband maser". Each single emission by a charged particle is due to the stimulated, 
as opposed to spontaneous, emission process (hereby the term maser). Unlike the conventional lasers in 
which basically one single frequency gets amplified, in this case charged particles can resonate with many 
mutually incoherent waves with different frequencies. 



In this paper we discuss a theory of pulsar radio emission developed by Lominadze, Machabeli 



k Mikhailovskii 1979, Machabeh fc Usov 1989, Kazbegi et al. 1991b, Kazbegi et al. 1991c, Lyutikov 



1998b, Lyutikov, Machabeli & Blandford 1998. We hypothesize, that pulsar radiation is generated by 
the instabilities developing in the outflowing plasma on the open field lines in the outer regions of the 
pulsar magnetosphere. Radiation is generated by two kinds of electromagnetic plasma instabilities - 
cyclotron-Cherenkov and Cherenkov-drift instabilities. The cyclotron-Cherenkov instability is responsible 
for the generation of the core-type emission and the Cherenkov-drift instability is responsible for the 
generation of the cone- type emission ( [Rankin 1990 ). The waves generated by these instabilities are 
vacuum-like electromagnetic waves so that they may leave magnetosphere the directly. 

In contrast to most modern theories of pulsar radio emission, cyclotron-Cherenkov and Cherenkov-drift 
instabilities occur in the outer parts of magnetosphere. The location of the emission region is determined 
by the corresponding resonant condition for the cyclotron-Cherenkov and Cherenkov-drift instabilities. 
Instabilities develop in a limited region on the open field lines. The size of the emission region is determined 
by the curvature of the magnetic field lines, which limits the length of the resonant wave-particle interaction. 
The location of the cyclotron-Cherenkov instability is restricted to those field lines with large radius of 
curvature, while the Cherenkov-drift instability occurs on field lines with curvature bounded both from 
above and from below. Thus, both instabilities produce narrow pulses, though they operate at radii where 
the opening angle of the open field lines is large. 

To a large extent a possible mechanism for the generation of pulsar radio emission is predicated on 
the choice of parameters of the plasma fiow that is generated by a rotating neutron star. At this point 
we know only the general features of the distribution function of the particles in a pulsar magnetosphere 
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(Tademaru 1973, Arons 1981b, Daugherty fc Harding 198^ ). It is believed to comprise (see Fig. |]) (i) a 



highly relativistic primary beam with the Lorentz factor 7b ~ 10^ and density equal to the Goldreich-Juhan 
density ngj = Jl • B/(2 tt e c), (ii) a secondary electron positron plasma with a bulk streaming Lorentz 
factor 7j, « 10 — 1000, a similar scatter in energy Tp « jp and a density much larger than the beam density 
Up w Xncj = 10'^ — lO^ficj, (iii) a tail of plasma distribution with the energy up to jt = 10* — 10^. 

The choice of a particular distribution function (Fig. |l|) is very important. All the following is 
dependent on this choice. In particular, in our model, the primary beam is composed of electrons or 
positrons. We show in Appendix ^ that both cyclotron-Cherenkov and Cherenkov-drift instabilities do not 
develop in an ion beam. 

The electromagnetic cyclotron-Cherenkov and Cherenkov-drift instabilities are the strongest instabilities 



in the pulsar magnetosphere ( Lyutikov 1997c ). This differs from the more common case of a nonrelativistic 
plasma, where electrostatic Cherenkov-type instabilities (i.e. those that result in emission of electrostatic 
Langmuir-type waves) are generally stronger than electromagnetic instabilities. In addition, for a 
one-dimensional plasma streaming along the magnetic field, the effective parallel mass is considerably 
increased by relativistic effects. For the particles in the primary beam, which contribute to the development 
of the instability, the effective parallel mass is rn^ffw — 10^^ to (to is a mass of a particle). This 

suppresses the development of the electrostatic instabilities. In contrast, the effective transverse mass , 
TOe//j_ = 7;,TO, is less affected by the large parallel momentum. The electromagnetic instabilities are less 
suppressed by the large streaming momenta. Thus, the relativistic velocities and one-dimensionality of 
the distribution function result in a strong suppression of the electrostatic instabilities and strengthen 
electromagnetic instabilities. 

Cyclotron-Cherenkov generation of wave by fast particles is not new in astrophysics. For example, 
cosmic rays in the interstellar medium and in supernova shocks generate Alfven by a similar mechanism. 
In the case of Alfven waves in the nonrelativistic electron-ion plasma, the frequency of the waves u) can 
be much smaller than the kv term and can be neglected in the resonance condition. The important 
difference between these applications and cyclotron-Cherenkov instability in pulsar magnetosphere is that 
the generated waves belong not to the hydromagnetic Alfven waves, that cannot leave the plasma, but to 
near vacuum electromagnetic waves. 

We should also mention that a cyclotron-Cherenkov instability of an electron beam propagating along 
a magnetic field is known in the laboratory as a very effective source of the high frequency microwave radio 



emission ( Galuzo et al. 1982 , Didenko et al. 1983 , Nusinovich et al. 1995 ). The so called slow- wave 



electron cyclotron masers (ECM) provides a high efhciency and high power microwave source. Though 
there are no commercial slow-wave ECM available now, they are believed to be very promising devices 
due to their better control of the beam quality and potentially more compact design than the cyclotron 
autoresonance masers. Thus, pulsars can be regarded as cosmic slow-wave ECMs. 



In previous work ( Lyutikov, Machabeli fc Blandford 1998 ) we developed a new approach to the 



amplification of curvature radiation. We argue, that a new, Cherenkov-drift instability, may be operational 
in the pulsar magnetosphere. The Cherenkov-drift emission combines the features of both Cherenkov and 
curvature emission processes. This instability is similar to the drift instabilities of the inhomogeneous 
plasma. The striking feature is that, unlike the nonrelativistic laboratory plasma, where drift instabilities 
develop on the low frequency waves with the wave length on the order of the inhomogeneity size, in the 
strong relativistic plasma drift instabilities can produce high frequency waves. We develop an approach that 
treat Cherenkov and curvature emission consistently in cylindrical coordinates. The choice of cylindrical 
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coordinatcs allows one to consider curvature emission as a resonant emission process, (in the former 
approaches the wave-particles interaction length was very limited, that precluded a strong amplification 
under all circumstances). Another important difference in this work is the proper account of the dispersion 
and polarization of the normal modes. We show that Cherenkov-drift instability develops only in a medium 
which supports subluminous waves. 

We argue that the theory based on the cyclotron-Cherenkov and Cherenkov-drift instabilities is 
capable of explaining a very broad range of the pulsars' observed properties. In a "standard" pulsar with a 
surface magnetic field B = 10^^ G and a period P = 0.5 s both cyclotron-Cherenkov and Cherenkov-drift 
instabilities occur at a radius of about 10^ cm. In a dipole geometry the emission region is limited to 
the field lines near the direction of the magnetic moment of the neutron star. The emission region for 
the Cherenkov-drift instability is larger than for the cyclotron-Cherenkov instability. In both cases it is 
determined by the curvature of the magnetic field lines that limit the coherent growth of the waves. 

In Section |^ we shortly review the observational properties of the pulsar radio emission. In Section 
^ we discuss the microphysics of the underlying cyclotron-Cherenkov and Cherenkov-drift instabilities. In 
Section ^ we describe the fiducial pulsar model, and finally in Section || we show how the properties of the 
pulsar radio emission may be explained in the framework of this theory. 



2. Observational properties and phenomenological theory of pulsar radio emission 

A useful observational framework for discussing theory is a description of pulsar radio emission given 



by Rankin 1990. The main feature of this model is the division of emission into two main classes: core and 
cone. There may be many cones of emission. In each pulsar the averaged profile may be a combination of 
core and/or cone emissions (Fig. |^). 

The majority of pulsars (about 70%) show core-type emission. The typical core emission has the 
following features: (a) the profile has a single component, (b) variable circular polarization (up to 60%), 
the amount of the circular polarization either reaches maximum at the maximum intensity or the sense of 
polarization may reverse in the middle of the pulse (Fig. emissiongeometry), (c) linear polarization changes 
from nearly 100% to unpolarized, in most cases the radiation may be split into two orthogonally polarized 



modes ( Stinebring et al. 1984) 



About 30% of all pulsars show pure conal emission and they are divided into two main groups - 
cone singles and cone doubles which are believed to be closely connected, the only difference being the 
geometrical path of the line of sight through emission region. Conal type emission shows a great variety 
of phenomena. Some of the typical features of conal emission are (a) the profiles can have up to four 
components, corresponding to two cones, (b) circular polarization is small and unsystematic, (c) linear 
polarization is moderate, some pulsars show a sudden change of position angle by ±7r/2; in most cases 
the radiation may be split into two orthogonally polarized modes, then the change of position angle by 
±7r/2 corresponds to the change in the relative intensity of the modes. Besides these phenomena, cone- type 
emission shows drifting subpulses, nuUings and mode switching. These effects are probably related to the 
temporal and/or spatial modulations of the parameters of the outflowing plasma and will not be discused 



here (see Kazbegi et al. 1996) 



There have been contradictory attempts to determine observationally the emission altitude but this 



conclusion is model-dependent. By contrast, Gwinn et al.l997 used interstellar scattering to measure 
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directly the size of the emission region of « 500 km QAc/fl). The other, model-dependent, observational 
evidences for high emission altitudes comes from the large duty cycles often observed in pulsars ("wide 
beam pulsars" 



e.g. 



Lyne fc Manchester 1 



Conventionally, these pulsars are interpreted as almost 



aligned rotators. An alternative explanation is that the emission may be coming from large radii. 

Another class of pulsars, those with interpulses, are conventionally interpreted as orthogonal rotators 
with emission coming from two poles despite the fact that emission bridges are often observed in these 



pulsars [Hankins fc Fowler 1986 and, at least in one case, PSR B0950-I-08 the polarization data imply a 
nearly aligned rotator ( Manchester 1995 ). The Crab pulsar also has a bridge of emission between the two 
pulses, intensities of the main pulse and interpulse are correlated and the angle of linear polarization in the 
two pulses seems to be related. In addition, it shows additional emission peaks between the interpulse and 



pulse at centimeter wave length Moffett & Hankins 1996. Manchester 1996 suggests that interpulses come 
from the same pole. If so, the simplest interpretation is that emission originates at high altitudes. 



3. Cyclotron-Cherenkov and Cherenkov-drift instabilities 



In this section we consider the physics of the cyclotron-Cherenkov and Cherenkov-drift instabilities 



( Ginzburg fc Eidman 1959 , Lyutikov, Machabeli fc Blandford 1998 ). The terminology used here to describe 
these instabilities refers to the fact that in the cyclotron-Cherenkov emission, a resonant particle changes its 
gyrational state (undergoes a transition between different Landau level) , thereby comes the " cyclotron" part 
of its name, but the force that induces the emission is due to the presence of a medium (the "Cherenkov" 
part of the name) . The Cherenkov-drift emission is similar to conventional Cherenkov (the gyrational state 
of the resonant particle remains unchanged) but it involves a nonvanishing curvature drift of the resonant 
particles. 

The interplay between cyclotron (or synchrotron) and Cherenkov radiation has been a long-standing 
matter of interest. 



Schwinger et al 1976 discussed the relation between these two seemingly different 



emission mechanisms. They showed that conventional synchrotron emission and Cherenkov radiation may 
be regarded as respectively limiting cases of |ri — 1| <C 1 and B = of a syncrgctic (using the terminology 



of Schwinger et al. Schwinger et al 1976 ) cyclotron-Cherenkov radiation. In the work ( Lyutikov, Machabeli 
fc Blandford 1998 this analogy has been extended to include inhomogeneity of the magnetic field. 



The physical origin of the emission in the case of Cherenkov-type and synchrotron-type processes is 
quite different. In the case of Cherenkov-type process, the emission may be attributed to the electromagnetic 
polarization shock front that develops in a dielectric medium due to the passage of a charged particle with 
speed larger than phase speed of waves in a medium. It is virtually a collective emission process. In the 
case of synchrotron-type process, the emission may be attributed to the Lorentz force acting on a particle 
in a magnetic field. Cherenkov-type emission is impossible in vacuum and in a medium with the refractive 
index smaller than unity. 

Both the cyclotron-Cherenkov and Cherenkov-drift instabilities, that we believe can develop in pulsar 



magnetosphere, operate in the kinetic regime, i.e. they are of a maser type (Lyutikov 1998a). This means 
that there is some kind of population inversion in the phase space, which supplies the energy for the 
development of the instability. In the present case, the source of free energy is the anisotropic distribution 
function of the fast particles. The condition of a population inversion may be restated that the induced 
emission dominates over induced absorption for a given transition. 
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The first two steps in identifying the possible maser-type radio emission generation mechanism are 
(i) determining which radiative transitions are allowed in a given system and (ii) establishing if the given 
distribution function allows for the population inversion for the particle in resonance with the emitted waves. 
In this chapter we will first discuss the microphysics of the two suggested emission mechanisms and then 
show that the distribution function of the particles present on the open field lines of pulsar magnetosphere 
does have a population inversion and allows maser action. When discussing the microphysics of the emission 
process, we will concentrate on the spontaneous emission processes. The induced emission rate, which is 
important for the development of the instabilities, is derivable from the spontaneous emission in the usual 
manner. In the process of induced emission, the electron emits a wave in phase with the incident wave. 
However, both cyclotron-Cherenkov and Cherenkov-drift masers are broadband and incoherent because a 
single electron can resonate with several waves simultaneously. 



3.1. Physics of Cyclotron-Cherenkov Emission 

The cyclotron-Cherenkov instability develops at the anomalous cyclotron resonance 

to(k) - k»v» + — = (1) 
7 

where Ci;(k) is the frequency of the normal mode, k is a wave vector, v is the velocity of the resonant 
particle, lub = \e\B/mc is the nonrelativistic gyrofrequency, 7 is the Lorentz factor in the pulsar frame, e 
is the charge of the resonant particle, m is its mass and c is the speed of light. Note a sign before the ujb 
term. 

To describe the microphysics of the cyclotron-Cherenkov emission emission process wc first recall the 
microphysics of the conventional Cherenkov emission ( pinzburg fc Eidman 1959 ). Consider a charged 



particle propagating in an unmagnetized dielectric with the dielectric constant e > 1. As the particle 
propagates, it induces a polarization in a medium. If the velocity of the particle is larger than the velocity 
of propagation of the polarization disturbances in a medium, which is equal to the phase speed of the waves 
Vph = c/^/e < c, the induced polarization cannot keep up with the particle. This results in a formation of 
the polarization shock front. At large distances, the electromagnetic fields from this "shock front" have a 
wave-like form corresponding to Cherenkov emission. Thus the emission is attributed to the polarization 
shock front and not directly to the particle. This polarization shock front acts on the particle with a drag 
force, which slows down the particle. This drag force may be considered as a generalization of the radiation 
reaction force in a medium. 

Now let us consider the propagation of a particle in a magnetized dielectric along a spiral trajectory. 
Similarly, the propagating charged particle induces polarization in a medium. If the velocity of the particle 
is larger than the phase speed of the waves, a polarization shock front develops, which acts on the particle 
with a drag force. Now the drag force, averaged over the gyrational period, has two components: along 
the external magnetic field and perpendicular to it. The parallel part of the drag force always slows the 
particle down. The surprising result is that the perpendicular component of the drag force acts to increase 
the transverse momentum of the particle. Thus a particle undergoes a transition to the state with higher 
transverse momentum and emits a photon. The energy is supplied by the parallel motion. 

The photons emitted by such mechanism correspond to the anomalous Doppler effect 
to — fc||W|| — slob/ Ires = 0, with s < 0. In a vacuum, only the normal Doppler resonance, with 
s > 0, is possible. The necessary condition for the anomalous Doppler effect, lo — fciiuii < 0, may be 
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satisfied for fast particles propagating in a medium with the refractive index larger than unity. It is natural 
to attribute the emission at the normal Doppler resonances to the Lorentz force of the magnetic field acting 
on the electron, while the emission at the anomalous Doppler resonances is attributed to the electromagnetic 
drag forces from the medium. 

It may also be instructive to consider the emission process in the frame associated with the center of 
gyration of the particle. In that frame the waves emitted at the anomalous Doppler effect have negative 
energy, so that in the emission process the particles increases its energy and emits a photon. 

The cyclotron-Cherenkov instability may be considered as a maser using the induced cyclotron- 
Cherenkov emission. The free energy for the growth of the instability comes from the nonequilibrium 
anisotropic distribution of fast particles. The condition that the emission rate dominates the absorption 
requires population inversion in the distribution function of fast particle (maser action). Since radiation 
reaction due to the anomalous Doppler effect induces transition up in quantum levels , for the instability to 
occur we need more particles on the lower levels. From the kinetic point of view, waves grow if the quantity 
k ^Q^^-* is positive for some values of k. For an electron in a magnetic field this condition takes the form 

sujb df , df 
v± dp± "dpii 

where s is a harmonic number. Here s > corresponds to normal Doppler effect (transition down in 
Landau levels) and s < corresponds to anomalous Doppler effect (transition up in Landau levels). If the 
distribution function is a plateau-like in parallel momenta then the condition for instability is s-^j- > 
which could be satisfied for inverted population for the normal Doppler effect or for the "direct" distribution 
for anomalous Doppler effect. The latter case takes place for the beam of particles propagating along the 
magnetic field with no dispersion in transverse momenta. 



3.2. Physics of Cherenkov-Drift Emission 

There is a possibility for the development of the Cherenkov-drift instability, which occurs at the 
resonance 

cj(k) — fc||ti|| — k±Ud = (3) 

where jv^^c/ lubRb is the drift velocity, Rb is the curvature radius of the magnetic field line. A weak 

inhomogeneity of the magnetic field results in a curvature drift motion of the particle perpendicular to 
the local plane of the magnetic field line. A gradient drift (proportional to (B • V)B much smaller than 
the curvature drift and will be neglected. When the motion of the particle parallel to the magnetic field 
is ultrarelativistic, the drift motion can become weakly relativistic even in a weakly inhomogeneous field 
resulting in the generation of electromagnetic, vacuumlike waves. The presence of three ingredients (strong 
but finite magnetic field, inhomogeneity of the field and a medium with the index of refraction larger than 
unity) is essential for this type of emission. We will call this mechanism Cherenkov-drift emission stressing 
the fact that microphysically it is virtually Cherenkov-type emission process. 



Conventional consideration of the curvature emission ( Blandford 1975, Zheleznyakov & Shaposhnikov 



197E , Luo fc Melrose 1992a , Melrose 1978b| ) emphasize the analogy between curvature emission and 



conventional cyclotron emission. To our opinion this approach, though formally correct, has limited 
applicability and misses some important physical properties of the emission mechanism. It has two 
important shortcomings. The first is that in adopting a plane wave formalism, the interaction length for an 
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individual electron, « Rs/jt, is essentially coextensive with the region over which the waves interact with 
a single electron. The approach necessarily precludes strong amplification under all circumstances because 
the wave would have to grow substantially during a single interaction in a manner that could not be easily 
quantified. The second problem was a neglect of dispersion. We address the first shortcoming by expanding 
the electromagnetic field in cylindrical waves centered on ?' = 0, and the second explicitly by considering 
general plasma modes. 



In a separate approach Kazbegi et al. [Kazbegi et al. 1991b| considered this process calculating a 



dielectric tensor of an inhomogeneous magnetized medium, thus treating the emission process as a collective 
effect. They showed that maser action is possible only if a medium supports subluminous waves. In 



the previous work (Lyutikov, Machabeli & Blandford 199S we showed how these two approaches can be 



reconciled and argue that the dielectric tensor approach, which treats the Cherenkov-drift emission as a 
collective process, has a wider applicability. 

It is more natural to consider Cherenkov-drift emission in a curved magnetic field as an analog of 
the Cherenkov emission with the drift of the resonance particles taken into account, than as the type of 
a curvature emission. From the microphysical point of view the emission is again due to the polarization 
shock front that develops due to the passage of a superluminal particle through a medium, so it is required 
that the emitting particle propagates with the velocity greater than the phase velocity of the emitted 
waves. The Cherenkov-drift maser is impossible in vacuum, unlike the curvature emission, which is a close 
analog of the conventional cyclotron emission and is possible in vacuum. The curvature provides only the 
drift component of the velocity, which is essential for the coupling the resonant particle to the emitted 
electromagnetic wave. 

In a Cherenkov-type emission the resonant particle can interact only with the part of the electric field 
parallel to the velocity. Thus, if the drift of the resonant particles perpendicular to the plane of the field 
line is taken into account, it becomes possible to emit a transverse electromagnetic wave with the electric 
field along the drift velocity, i.e. perpendicular to the plane of the curved field line (see Fig. ^). The 
growth occurs on the rising part of the parallel distribution function where > 0. This is satisfied for the 
particles of the primary beam. 



4. Model of the Pulsar Radio Emission 

4.1. Plasma generation 

Rotating, strongly magnetized neutron stars induce strong electric fields that pull the charges from 
their surfaces. Inside the closed field lines of the neutron star magnetosphere, a steady charge distribution 
established, compensating the induced electric field. On the open field lines, the neutron star generates a 
dense fiow of relativistic electron-positron pairs penetrated by a highly relativistic electron or positron beam. 
The density of the primary beam is roughly equal to the Goldreich- Julian density uqj = Jl • B/(2 tt e c). 
We will normalize the density of the pair plasma to the Goldreich-Julian density. 

Tip = Ancj = 10^ - lO^ncj, uj^ = Xlu^ = 2Xljb^ (4) 
where A is the multiplicity factor which is the number of pairs produced by each primary particle. Secondary 



pairs are born with almost the same energy in the avalanche-like process above in the polar cap (Arone 



1983). The pair creation front in the polar cap region is expected to be very thin so we can in the first 



approximation neglect the residual electric field in the front that could lead to the reversed current and 
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different initial energies and densities of secondary particles. The combination of the pair plasma and 
primary beam is expected to screen the rotationally induced electric field so that the flow is force-free. 

Another relation between the parameters of the plasma and the beam comes from the energy argument 
that the primary particles stop producing the pairs when the energy in the pair plasma becomes equal to 
the energy in the primary beam: 

2 < 7 >^'' A ~ 7f,, at the pair formation front (5) 

where it was assumed that the initial densities, temperatures and velocities of the plasma components are 
equal. For cold components < 7 >j?^ = 7^ while for the relativistic components with a temperature Tp the 
average energy is < 7 >'±' — '2,^pTp. The assumption of equipartition (|^) is very approximate one, but it 
allows considerable simplification of numerical estimates. If for some reason, this would turn out to be an 
incorrect assumption, the corresponding formula can be adjusted by changing the scaling. 

As an estimate of the densities of the particles from the tail of plasma distribution we will use the 
assumption that the energy in the tail is approximately equal to the energy in the plasma (and in the 
beam): 

Itnt « 7^°^ (6) 
where 74 and nt are the typical energy and the density of the tail particles. 



4.2. Typical pulsar 

In this work we will make numerical estimates for the "typical" pulsar with the following parameters: 

(i) magnetic field is assumed to be dipole with the magnetic field strength at the surface of neutron 
star Bns = lO^^G, 

(ii) rotational period of the star P = 0.5 s (light cylinder radius Rls = 2.4 x 10^ cm). 

(iii) average streaming energy of the plasma components 7^ = 10 

(iv) temperature of the plasma components Tp = 10. 

(v) initial energy of the primary beam at the pair formation front 7^ = 6 x 10^. 

The energy of the beam will decrease due to the curvature radiation reaction force (Appendix ^) . Then 
at the light cylinder, where the instabilities occur, the beam will decelerate due to the curvature radiation 
reaction to 7^ = 2 x 10^. 

For a given period and magnetic field the equation (^) reduces the number of free parameters for 
the plasma to two: the plasma temperature and the bulk streaming energy jp (or temperature and the 
multiplicity factor A). We chose a strongly relativistic plasma with the invariant temperature Tp = 10. The 
multiplicity factor A corresponding to these parameters follows from Eq. (||): A = 3 x 10^. The average 
energy of the tail particles is assumed to be 74 = 10^. An important factor that determines the growth 
rate of the instabilities is the energy scatter of the resonant particles. In estimating the growth rates of the 
cyclotron-Cherenkov and Cherenkov-drift instabilities on the primary beam we will also assume that the 
scatter in Lorentz factors of the primary particles in the pulsar frame A7 = 10^. This assumes that the 
beam has cooled considerably due to the curvature radiation and lost about 90% percent of its energy. 



- 10 - 



Several points are important in our choice of parameters. First, we use a relatively low plasma 



streaming 7-factor (and respectively high multiplicity factor A). In the polar cap models ( Arons 1981a 



Arons 1982) the pulsar plasma will have a low streaming 7-factor if the the magnetic field near the surface 



differs considerably from the dipole field thus reducing the radius of curvature ( {Machabeli fc Usov 1989 ) 



Secondly, the required scatter in energy of the primary beam particles (A7 — 10^) is very small. This is 
due to the effects of curvature radiation reaction on the primary beam during its propagation through the 
dipole pulsar magnetosphere. The highly nonlinear damping rate due the emission of curvature radiation 
by the primary beam result in an effective cooling of the beam (see Appendix^). 

We used two types of the particle distribution function to calculate the relevant moments: water bag 



and relativistic Maxwellian (Lyutikov 1998a). For the case of relativistic Maxwellian distribution function 
the strongly relativistic temperature Tp — 10 implies that, formally, there are many backward streaming 
particles. We note, though, that the backward streaming particles do not contribute significantly to any of 
the relevant moments of the distribution function, so that we can regard the strongly relativistic streaming 
Maxwellian distribution as a fair approximation to the relatively unknown, but definitely very hot, real 
distribution function. 

The location of our typical pulsar on the P — E diagram is shown on Fig ^. The corresponding plasma 
densities and frequencies are given in Tables |l] and || for the two locations (near the stellar surface and at 
the emission region R sa 10^ cm in the pulsar and plasma reference frames. 

The radial dependence of the parameters is assumed to follow the dipole geometry of the magnetic 

field: 

UJBir) = u;b{Rns) f 

ujp{r) ^ ujp{RNs) {—^\ ■ (7) 

This may not be a good approximation in the outer regions of the pulsar magnetosphere, where relativistic 
retardation, currents flowing in the magnetosphere and the effects of plasma loading considerably change 
the structure of the magnetic field. 



4.3. Fundamental plasma modes 

It is important for the development of both instabilities that the plasma supported subluminous waves. 
Such wave indeed exist in a strongly magnetize electron-positron plasma. If we neglect (as a first order 
approximation) the relative streaming of the plasma electrons and positrons and the field curvature then 
there are three fundamental modes: transverse extraordinary wave with the electric field perpendicular 
to the k-B plane and longitudinal-transverse mode the electric field in the k-B with two branches: 



ordinary mode and Alfvcn mode (Arons & Barnard 1986, Volokitin, Krasnosel'skikh & Machabeli 1985 



Lyutikov 1998a). For small angles of propagation with respect to magnetic field the ordinary mode (upper 
longitudinal-transverse wave) is quasi-longitudinal for small wave vector (fc <C Wp/c, Up - plasma frequency) 
and quasi-transverse for large wave vector (fc ^ Wp/c), while the Alfven mode is quasi-transverse for small 
wave vector and quasi-longitudinal for large ones. Fig. ^ illustrates these modes in the case of cold plasma 
components in the plasma frame and in the low frequency limit uj' ujb (w' is the wave frequency in the 
plasma frame). 





Pulsar frame 


Plasma frame 


iVldi^lictiC ilcitlj Vjr 






Cyclotron frequency Wb*, (rads~i), 


1.8 X 10^'' 


1.8 X lOi'' 


Beam density, cm~^ 


n'cj = 05/(2 Tree) = 1.4 X lO^i 


no J = f2B/(27rec7p) = 1.4 x 10i° 


Beam plasma frequency, rads~i 


^'gj = ^'-^^ = V20a;B = 2 x lO^o 


ujGj = ^ = = 6.3 x 109 


Beam energy 


7^ = 6 X 10^ 


76 = 2^ = 3 X 106 


Plasma density, cm~^ 


n' = ^''^ = 4 X 10^6 

P 2 TT e c 


"P'= 2^=4x1015 


Plasma frequency, rads~i 


= V2wbOA = 1.2 X 10^^ 


cop = ^^-'Bnx ^ 4 X 1012 


1 — Px, Px- phase speed of X mode 


S' = ff^ - 2^ = 10-15 





Table 1: Plasma parameters at the surface of the neutron star 





Pulsar frame 


Plasma frame 


Magnetic field, G 


10^ 


10^ 


Cyclotron frequency * , (rads~^), 


1.8 X 10^° 


1.8 X 10^° 


Beam density, cm~^ 


n'dj = nB/(2 7rec) = 1.4 x 10^ 


no J = / {2 TT e cjp) = 14 


Beam plasma frequency, rads~^ 


u;'^j = = 6X 105 


WGJ = # = 3 X 105 


Beam energy 


7^ « xlO^ 




Plasma density, cm~^ 


"p 27rec ^ X -^^ 


«P'= 2^=4x106 


Plasma frequency, rads~^ 


U'p = V2WBiU = 4 X 10° 


UJp = ^2a,BnA ^ ^ 2 X 10^ 


1 — Px,f3x phase speed of X mode 


6' = = = 5 X 10-^ 


S = ^ = ^^"^p = 1.2 X 10-^ 


Typical frequency, rads~^ 


w' = 5 X 10^ 


w = ^ = 2.5 X 10^ 



Table 2: Plasma parameters at the emission region ii 10® cm 
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The dispersion relations for the extraordinary mode in the pulsar frame in the limit uj <ti is 

L.^kc{l-S), S^^^ (8) 

where we used a relationship between the plasma density and the Goldrcich- Julian density Tip — Xtiqj and 
the definition of the Goldreich-Julian density. 

For the coupled Alfven and ordinary modes it is possible to obtain the asymptotic expansion of the 
dispersion relations in the limits of very small and very large wave vectors. Here we will give the dispersion 
relations for the quasitransverse parts of the waves 



(l - + '^"""5""'" ) 0-wave if kc » ^pipUjp and 9 < 1/Tp 

c2 fc2 cos2 0(1- -4^ - ,"'^'""'0 Alfven wave if kc < J%-lr,LOr, (9) 



These relationships are valid for the frequencies satisfying the inequality 

uj < -ipUs/Tp (10) 

This is a condition that in the reference frame of the plasma the frequency of the waves is much smaller 
that the typical cyclotron frequency of the particles cj^/Tp. 

When the relative streaming of the plasma is taken into account, the dispersion relations change 
considerably: a new, slow Alfven branch appears and the extraordinary and Alfven modes behave differently 



in the region of small wave vectors (Lyutikov 1998b). The polarization of the fundamental modes also 



changes. For the angles of propagation with respect to magnetic field less than some critical angle, 
which depends on the difference of the velocities of the components and on the magnetic field, the two 
quasi-transverse waves are circular-polarized, while for the larger angles the two transverse modes become 
linearly-polarized. 



4.4. Development of instabilities 

In this section it is shown, that the pulsar radiation may be generated by two kinds of electromagnetic 
plasma instabilities - cyclotron-Cherenkov and Cherenkov-drift instabilities. The cyclotron-Cherenkov 
instability is responsible for the generation of the core-type emission and the Cherenkov-drift instability 



is responsible for the generation of the cone-type emission (Rankin 1990). The waves generated by these 
instabilities are vacuum-like electromagnetic waves: they may leave the magnetosphere directly. 

We assume that the distribution function, displayed in Fig. |l|, remains unchanged throughout the 
inner magnetosphere. This requires that no Cherenkov-type two-stream instabilities develop and that the 
high energy particles are not excited to the high gyrational states by the mutual collisions or by the inverse 
Compton effect. Then the outer regions of the pulsar magnetosphere two instabilities can develop: (i) 
cyclotron-Cherenkov instability and (ii) Cherenkov-drift instability. 

A detailed consideration of the conditions necessary for the development of the cyclotron-Cherenkov and 
Cherenkov-drift instabilities are given in Appendix ^ and |C[ Both cyclotron-Cherenkov and Cherenkov-drift 
instabilities develop in the outer regions of the pulsar magnetosphere at radii R ~ 10^ cm. 
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The frequencies of the waves generated by the cyclotron-Cherenkov instabihty are given by ( [A3| ) 

^1 



473^1 



Ti 2 

which may be solved for the radius at which the waves with frequency oj are emitted 



(11) 



R-Rns[ ] f ) = 1.8 X 10^ Bli' .-/^ Po^iS-l^^m (12) 



The relationship (12) may be regarded as a " radius-to- frequency" mapping. For a given jres the radial 
dependence of the right hand side of equation (M ) will result in a radial dependence of the emitted frequency. 
The radial dependence of the parameters in ( [A3 ) will result in emission of higher frequencies deeper in the 
pulsar magnetosphere, exactly what is observed. The frequencies emitted at the Cherenkov-drift resonance 
do not have a simple dependence on radius from the neutron star. They are determined by several emission 
conditions which limit the development of the instabilities to the particular location in the magnetosphere 
and particular frequencies. 

The conditions that the instabilities should satisfy are: 

(i) small growth length in a curved field lines c/F < A9Rb, where A6 is the range of the emitted 
resonant angles 

(ii) condition of kinetic instability |k(5vres| S> F 

In addition to these, the Cherenkov-drift instability should also satisfy another condition: 

(iii) the condition of a large drift Ud/c> V^d. 

The condition (i) states that an emitting particle can stay in a resonances with the wave for many 
growth lengths. The condition (ii) is a requirement that the growth rate of instability is much smaller than 
the bandwidth of the growing waves. This condition is necessary for the random phase approximation to 
the wave-particle interaction to apply. 

In Appendices |a[ |c| and ^ we show that the above conditions can be satisfied for the chosen set of 
parameters both in normal pulsar and millisecond pulsars. The conditions for the development of the 
cyclotron-Cherenkov and Cherenkov-drift instabilities depend in a different ways on the parameters of the 
plasma. They may develop in the different regions of the pulsar magnetosphere. 



5. Pulsar Phenomenology 
5.1. Energetics 

In both cyclotron-Cherenkov and Cherenkov-drift mechanism, the emission is generated by the fast 
particles which supply the energy for the growth of the waves. The total energy available for the conversion 
into radio emission is of the order of the energy of the particle flow along the open field lines of pulsar 
magnetosphere: 

« nGj7ri?2^7,,mc3 « lO^^ergs-i (13) 

where Rpc — RNS^J^^^ is the radius of a polar cap. This is sufficient to explain the radio luminosities 
of the typical pulsar if the effective emitting area is about one hundredth of the total open field line cross 
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section. 



5.2. Emission Pattern 



The emission pattern for the "core" -type pulsars (according to the classification of (Rankin 199C) is 
a circle with the angular extent of several degrees. In our model the region of the cyclotron-Cherenkov 
instability is limited to the nearly straight field lines (see Fig. ^). 

The curvature of the magnetic field destroys the coherence between the waves and the resonant 
particles. To produce an observable emission the waves need to travel in resonance with the particle at 
least several instability growth lengths. Far from the magnetic axis, where the curvature is substantial, the 
waves leave the resonance cone before they travel a growth length and no substantial amplification occurs. 
Near the magnetic field axis the radius of curvature is very large and waves can stay in a resonance with a 
particles for a long time and grow to large amplitudes. 

By contrast, the Cherenkov-drift instability requires curvature of the field lines, but its growth rate may 
be limited by the coherence condition. In a dipole geometry, this will limit the Cherenkov-drift emission 
region to a ring around the central field line. Another possible location of the Cherenkov-drift instability is 
the region of the swept back field lines (Fig. |^ 

The development of the cyclotron-Cherenkov instability depends on fewer parameters than the 
Cherenkov-drift instability. It develops on the nearly straight central field lines. The conditions for the 
development of the cyclotron-Cherenkov instability (Eqs. A4, A9, |A13| ) depend on low powers of the 



plasma parameters, so it is quite robust. By contrast, the Cherenkov-drift instability depends on the 
plasma parameters and the radius of curvature of the magnetic field in a complicated way. This results in a 
broader range of phenomena observed in the cone emission. If the parameters of the plasma change due the 
changing conditions at the pair production front, the location of the Cherenkov-drift instability may change 
considerably. This may account for the mode switching observed in the cone emission of some pulsars. 

It may also possible beto explain in the framework of our model the phenomenon of "wide beam 



geometry" observed in some pulsars (Manchester 1996). The Cherenkov-drift instability may occur in the 



region, where the field lines are swept back considerably. Then the emission will be generated in what could 
be called a "wide beam geometry". 



5.3. Polarization 



If the average energy of the electrons and positrons of the secondary plasma is the same, the fundamental 
modes of such strongly magnetized plasma are linearly polarized. Both of the two quasi-transverse modes 
(one with electric field lying in the k — Bq plane, another with E perpendicular to this plan) may be 
emitted by the cyclotron-Cherenkov mechanism. This may naturally explain the two orthogonal modes 
observed in pulsars (Kazbegi et al. 1991b, Kazbegi et al. 1991c, Lominadze, Machabeli & Usov 1983). The 
difference in the dispersion relations and in the emission and absorption conditions between ordinary and 
extraordinary modes may explain the difference in the observed intensities of these modes. 



The rotation of the magnetized neutron star produces a difference in the average streaming velocities 
of the plasma components. This results in the circular polarization of the quasi-transverse modes for the 
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angles of propagation with respect to the local magnetic field less than some critical angle. This may 
explain in a natural way the occurrence of the circular polarization in the " core" emission. 

If the cyclotron-Cherenkov instability occurs on the tail of the plasma distribution function then the 
particles of both signs of charge can resonate with the wave. In a curved magnetic field, electrons and 
positrons will drift in opposite directions (Fig. ^). As the line of sight crosses the emission region, the 
observer will first see the left circularly polarized wave emitted by the electrons in the direction of their 
drift. When the line of sight becomes parallel to the local magnetic field the wave will resonate with both 
electrons and positrons in the plasma tail, so that the resulting circular polarization will be zero. Finally, 
the observer will see the wave emitted by positrons in the direction of their drift. This can explain the 
switch of the circular polarization observed in some pulsars. 

The cone emission, which is due to the Cherenkov-drift instability naturally has one linear polarization. 
An important difference from the standard bunching theory is that the waves emitted at the Cherenkov-drift 
resonance are polarized perpendicular to the plane of the curved magnetic field line. This may be used 
as a test to distinguish between the two theories. To do so, one need to determine the absolute position 
of the rotation axis of a pulsar - a possible but a difficult task. One possible experiment would involve 



Harrison- Tademaru effect (Harrison & Tademaru 1975), which predicts that the spins of neutron stars may 



be aligned with their proper motion due to the quadrupole magnetic radiation if the magnetic moment is 
displaced from the center of the star. Unfortunately, the current data does not support this theory (in the 
Crab, the spin of the neutron star is within 10° of from the direction of the proper motion, while for Vela it 
is approximately perpendicular). Symmetry of plerions and direct observation of jets from pulsars (like the 
one observed in Vela pulsar ( Markwardt fc Ogelman 1995| )) combined with polarimetry of the pulsar may 



be useful in determining the relative position of the electric field of the wave and the magnetic axis. 



5.4. Radius-to-frequency mapping 

The cyclotron-Cherenkov mechanism predicts a simple dependence of the emitted frequencies on 
the altitude (Eq. ^2|). The predicted power-law scaling of the emission altitudes is i?(i') oc This 



is strikingly close to the observed scaling cx v i^^ i ( Lesch et al. 199S , Kramer et al. 1994 ), which is 



derived from the simultaneous multifrequency observations of the time of arrival of pulses. A simple 

" radius-to- frequency" mapping will be "blurred" by the scatter in energies of the resonant particles, but the 

general trend that lower frequencies are emitted higher in the magnetosphere will remain. 

The Cherenkov-drift instability, on the other hand, does not have a simple dependence of the emission 
altitudes on the frequency. The resonance conditions for the Cherenkov-drift instability are virtually 
independent of frequency (Eq. pi] ), so that the location of the emission region is determined by the various 
conditions on the development of the instability (Appendix 



5.5. Formation of Spectra 

Development of the Cyclotron-Cherenkov and Cherenkov-drift instabilities results in an exponential 
growth of the electromagnetic waves. The original growth may be limited by several factors. A very likely 
possibility is that a spatial growth of the waves is limited by the changing parameters of the plasma. This 
possibility is hard to quantify since we do not know structure of the magnetic field in the outer regions of 
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the pulsar magnetosphere. The original growth of the instabilities may also be limited by the nonlinear 
processes: quasi-linear diffusion, induced scattering and wave decay. Finally, as the waves propagate in a 



magnetosphere, they may be absorbed by the particles of the bulk plasma (Lyutikov 1998a). The emergent 



spectra are the combined results of the these processes: emission, nonlinear saturation and absorption. 
We have considered the two most likely nonlinear saturation effects for the cyclotron-Cherenkov 



instability: quasilinear diffusion and induced Raman scattering (Lyutikov 1998c, Lyutikov 1998d). In 
the case of quasilinear diffusion, the induced transitions to upper Landau levels due the development 
of the cyclotron-Cherenkov instability is balanced by the radiation reaction force due to the cyclotron 
emission at the normal Doppler resonance and the force arising in the inhomogeneous magnetic field due 
to the conservation of the adiabatic invariant. These forces result in a quasilinear state and a saturation 
of the quasilinear diffusion. We have found a state, in which the particles are constantly slowing down 
their parallel motion, mainly due to the component along magnetic field of the radiation reaction force of 
emission at the anomalous Doppler resonance. At the same time they keep the pitch angle almost constant 
due to the balance of the force arising in the inhomogeneous magnetic field due to the conservation of the 
adiabatic invariant and the component perpendicular to the magnetic field of the radiation reaction force 
of emission at the anomalous Doppler resonance. We calculated the distribution function and the wave 



intensities for such quasilinear state (Lyutikov 1998d) 



In the process of the quasilinear diffusion, the initial beam looses a large fraction of its initial energy 
« 10%, which is enough to explain the typical luminosities of pulsars. The theory predicts a spectral index 
F{v) cx is the spectral flux density) which is very close to the observed mean spectral index of 



— 1.6 (Lorimcr at al. 1995). The predicted spectra also show a turn off at the low frequencies v < 300MHz 
and a flattering of the spectrum at large frequencies ly > IGHz (this may be related to the upturn in pulsar 
spectra observed at mm- wavelength, ( Kramer et al. 1996| )). 



The other possible mechanism that may be important for wave propagation and as an effective saturation 
mechanism for instabilities of electromagnetic waves is induced Raman scattering of electromagnetic waves 



(Lyutikov 1998c). The frequencies, at which strong Raman scattering occur in the outer parts of 
magnetosphere, fall into the observed radio band. The typical threshold intensities for the strong Raman 
scattering are of the order of the observed intensities, implying that pulsar magnetosphere may be optically 
thick to Raman scattering of electromagnetic waves. 



Absorption processes can play an important role in the formation of the emergent spectra (Lyutikov 



1998a). The waves may be strongly damped on the three possible resonances: Cherenkov, Cherenkov-drift 
and cyclotron. Alfven wave is always strongly damped on the Cherenkov resonances and possibly on the 
Cherenkov-drift resonance and cannot leave magnetosphere. Both ordinary and extraordinary waves may be 
damped on the Cherenkov-drift resonance. In this case Cherenkov-drift resonance affects only the ordinary 
and extraordinary waves with the electric field perpendicular to the plane of the curved magnetic field line. 
The high frequency parts of the extraordinary and ordinary wave may also be damped on the cyclotron 
resonances. In fact, in a dipole geometry, electromagnetic waves propagating outward are always absorbed. 
The fact that we do see radiation implies that waves get " detached" from plasma escaping absorption. 



5.6. High Energy Emission 



It may be possible also to relate the pulsed high energy emission to the radio mechanisms. The 
development of the cyclotron-Cherenkov instability at the anomalous Doppler effect leads to the finite pitch 



- 18 - 



angles of the resonant particles. The particles will undergo a cyclotron transition at the normal Doppler 
effect decreasing the their pitch angle. The frequency of the wave emitted in such a transition will fall in 
the soft X-ray range with the frequencies 

uj « 7f,tjB = lO^^rads"^ {E « IkeV) (14) 

In such a model the high energy emission will coincide with the core component of the radio emission, which 
is what is observed in Crab pulsar and some other pulsars. We also note, that this may be a feasible theory 
for the pulsed soft X-ray emission. The hard X-ray and 7-emission cannot be explained by this mechanism 
since the total energy flux in the primary beam is not enough to account for the very high energy emission 



(see, for example, Usov 1996) 



6. Observational Predictions 
General Predictions of the Maser-type Instability 

The formalism of the maser-type emission assumes a random phase approximation for the wave-particle 
interaction. Any given particle can simultaneously resonate with several waves whose phases are not 
correlated. This is different from the reactive-type emission (like coherent emission by bunches) when the 
intensities of the waves with different frequencies are strongly correlated. Thus any observation of the 
fine frequency structure in the pulsar radio emission may be considered as a strong argument against the 
reactive-type emission and in favor of the maser-type emission. 

Frequency Dependence of the Circular Polarization in the Core 

The cyclotron- Cherenkov instability can develop both on the primary beam particles and the particles 
from the tail of the plasma distribution. The cyclotron-Chcrenkov instability on the primary beam produces 
a pulse which has a maximum circular polarization in its center. The cyclotron- Cherenkov instability on 
the tail particles produces a pulse with the switching of the sense of the circular polarization in its center. 



Since the resonance on the tail particles occurs on larger frequencies ( Eq. A3 ) the effects of the switching 



sense of the circular polarization should be more prominent on the higher frequencies. 
Linear Polarization of the Cone Radiation 

Cherenkov-drift instability produces waves with the linear polarization perpendicular to the plane of 
the curved field line. This is in a sharp contrast to many other theories of the radio emission that tend to 
generate waves with the electric field in the plane of the curved field line. If one can determine the absolute 
position of the rotational axis, magnetic moment and the electric vector of the linear polarization, then, 
assuming a dipole geometry, it will be possible to determine unambiguously the position of the electric 
vector of the emitted wave with respect to the plane of the magnetic field line. 

No Cyclotron-Cherenkov Instability in Millisecond Pulsars 

In Appendix ^ we show that the cyclotron-Cherenkov instability does not develop in the magnetospheres 
of the millisecond pulsars. Since in our model the cyclotron-Cherenkov instability produces core- type 
emission we predict that millisecond pulsars will not show core- type emission. At present, the existing 
observations do not allow a separation into the core and cone type emission in the millisecond pulsars. 
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7. Conclusion 

In this work we presented a model for the pulsar emission generation capable of explaining a broad 
range of the; observations including the morphology, polarization and spectrum formation in normal and 
millisecond pulsars. We provided no explanation for the temporal behavior, nulling, drifting subpulses 
although these phenomena do not challenge the model. Our model makes several testable predictions, like 
high altitude of the emission region and an unusual relation of the polarization direction to the magnetic 
axis for the cone emission. In addition it has flexibility to accommodate a variety of phenomena arising, 
for example, due to the temporal variations in the flow and the structure of the magnetic field in a given 
pulsar, or due to the different structure of the magnetosphere in different pulsars. 

Further progress requires a better understanding of the structure of the pulsar magnetospheres both 
near the stellar surface where the mutipole moments of the magnetic field can substantially affect the 
physics of the pair formation front, and near the light cylinder, where the cmiission is taking place. It is also 
necessary to consider the nonlinear stages of the emission mechanisms described above. 

We would like to thank George Melikidze for his comments. ML would like to thank the Abastumani 
Astrophysics Observatory for the hospitality during his stays in Tbilisi and GZM acknowledges the support 
and hospitality during his stay at Caltech. This research was supported by NSF grant AST-9529170. 
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A. Conditions on cyclotron-Cherenkov instability 

In this Appendix we consider the development of the cyclotron-Cherenkov instability in magnetosphere 
of the typical pulsar. We show that the conditions on page |lj (Section 4.4) for the development of the 



cyclotron-Cherenkov instability instability are satisfied for the typical pulsar. 

The conditions for the development of the cyclotron-Cherenkov instability may be easily derived for 
the small angles of propagation with respect to the magnetic field. Representing the wave's dispersion as 

X mode 



47|^ 2 2 

u; = fcc(l - S), where S = { - '''^lll f O mode (Al) 

-4^ -t- 'f'^'^'f^'^ Alfvcn mode 

and neglecting the drift term, the resonance condition may then be written as 

2 "/res 2 LU -/res 

Let us discuss the condition for the development of the cyclotron-Cherenkov instability. First we note 



that equation (A2) requires that ^ \ < 5 and ^ < 6. The first is the condition that the particle is 
moving through plasma with the velocity faster than the phase velocity of the wave. The second condition 
limits the emission to the small angles with respect to magnetic field. Assuming that <C 5 and 



^ (5 we find from (A2) 



- = ^ (A3) 

IresO 



This resonant frequency increases with radius as . 

Using the upper limit on the frequencies for the relations (|Al[ ) to hold, namely oj <^ ^pLOB/Tp. This 
sets the limit on 5: 5 ^ Tp/ {^p^res)- This condition is satisfied for the radii: 

( R \ ( 2 7p^B,jvs \ ^'^ _ { 5 X 10^ for the beam 

\Rns ) ^ \ JbXn j ^ \ 2 X 10^ for the tail *■ ' 

This is the radius where the resonance condition first becomes satisfied for the particles of the primary 
beam and from the tail of the distribution function. 



The cyclotron-Cherenkov instability growth rate is (e.g., Kazbegi et al. 1991c): 

(A5) 



^p,res / X ^res^p ^res ^ 



2 LU A7 V 2 A7 7^^ LUb 

where we have normalized the density of the resonant particles to the Goldreich-Julian density 

, i2 — \ , ,2 
^p.res ~ '^res lUqj- 

It follows from ( [A5| ) that the growth rate increases with radius as R^. Deeper in the magnetosphere 
the growth rate is slow and the waves are not excited. At some point the growth rate becomes comparable 
to the dynamic time F/il = 1. This occurs at 

X 1/3 

(A6) 



f ^^^l ^*b\ _l 2x10^ for the beam 

\Rns) ^ \ ^\KesTpires 1x10^ for the tail 
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Using the equipartition condition we conclude that at a given radius the growth on the tail particles 
is approximately the same (the higher resonant density is compensated by the higher resonant 7-factor. 

Starting this radius the waves will start to grow with the growth rate increasing with radius. The 



highest frequency of the growing mode is determined by the condition ( [A3| ) evaluated at the radius (A6): 

1 X lO^rads^^ for the resonance on the beam /'A7^ 
4 X lO^^rads^^ for the resonance on the tail 

These estimates show, that for the chosen plasma parameters the cyclotron-Cherenkov instability always 
develops on the tail of the distribution function and can also develop on the particles from the primary 
beam. The higher density of the tail particles favors the development of the instability on the tail particles. 
The cyclotron-Cherenkov instability on the tail particles occurs deeper in the magnetosphere on the higher 
frequencies than the cyclotron-Cherenkov instability on the primary beam, which can develop further out 
in the magnetosphere and produce emission at the lower frequencies. 

The growth rate of the cyclotron-Cherenkov instability should satisfy two other conditions: the growth 
length must be much less that the length of the coherent wave-particle interaction and the condition of the 
kinetic approximation. The coherent wave-particle interaction in limited by the curvature of the magnetic 
field. A particle can resonate with the waves propagating in a limited range of angles with respect to the 
magnetic field. As the wave propagates in a curved field the angle that the wave vector makes with the 
field changes and the wave may leave the range of resonant angles. If the range of the resonant angles is 
then the condition of a short growth length is 

^ < AORb (A8) 

For cyclotron-Cherenkov instability we can estimate Ad « \/S. Then we find a condition on the radius 
of curvature 

1 cA'j ( lO^^cm for the beam /An^ 

^ ^ ^res^lresS^/^ \ 5 X lO^^cm for the tail ^ ' 

at the distance i? = 2 x lO^cm. The region of the cyclotron-Cherenkov instability is limited to the 
field lines near the central field line. The transverse size of the emission region may be estimated as 
X = R^/Rb ~ lO^cm. This gives an opening angle of the emission 0"^™ — xj Rk, 2°. 



There is another condition that the growth rate (A5) should satisfy, namely the condition of the kinetic 



approximation. In deriving the growth rate we implicitly assumed that the wave-particle interaction is 
described by the random phase approximation, which requires that the spread of the resonant particles 
satisfies the condition 

|k<5v,,,|»r (AlO) 
For the particle streaming in the curved magnetic field without any gyration this condition takes the 

form 



A7 A7 sujbA^ 



>r (All) 



For the cyclotron-Cherenkov instability the last term on the left hand side of (All) is the dominant: 

LOB A7 



^2 

Ires 



> r (A12) 
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Using the growth rate ( |A5| ) this condition may be rewritten as 

1/6 




3 X 10^ (A13) 



Equation ( A13| ) means that the kinetic approximation is weU satisfied inside the pulsar magnetosphere. 



Thus, the conditions for the development of the cyclotron-Cherenkov instability may be satisfied for 
the typical pulsar parameters. The waves are generated in the observed frequency range near the central 
field line. The radius of emission is ~ lO^cm, the transverse size of the emitting region is ~ lO^cm and the 
"thickness" of the emitting region is < lO^cm. This may account for the core- type emission pattern. 



B. Grov^r rate for the Cherenkov-drift instability 



In this Appendix we calculate the growth rate for the Cherenkov-drift instability in the plane wave 



approximation. Originally this growth rate has been obtained by Lyutikov, Machabeli & Blandford 199S 



using the single particle emissivities. Here we give a derivation using the simplified dielectric tensor 



(Lyutikov, Machabeli & Blandford 1998). The relevant components of the antihermitian part of the 
dielectric tensor are ( Kazbegi et al. 1991b ) . 



4^2^2 

i 

UJC 



dP4> 



OJC 



dfip^) 



5 (lu — I 



dp4. 



kxUd) , 
- k^Ud) 
kxUd) ■ 



(Bl) 



The growth rate of the Cherenkov drift instability can be obtained from (Melrose 1978b) 

a/(p) 



r = y dpu;(k, p) Kk ■ 



9p 



The growth rate for the It-mode is 



ckkj_ c k 



dprj, 



and growth rate for the t-mode is 



dp4, 



k±c J dp^ 



6{uj - k^v^ - kxUd) 



(5 (w - k^v^ - kxUd) ■ 



(B2) 



(B3) 



(B4) 



Where we chose the following polarization vectors 

1 



6 — kk i^'^^'' kj_^kxkfp^ 
e =- {—kx,0,kr,}, 

fc_L 



(B5) 
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where kj_ = yjk'^ + k'^ and 

k = {fcr,/c0,/c3.}. (B6) 



(see Lyutikov, Machabeli fc Blandford 1998 ) for the discussion of the polarization properties of normal 
modes in anisotropic dielectric in cylindrical coordinates. 

The maximum growth rate for the t-mode is reached when kx/k^ = Ud/c and the maximum growth 
rate for the It-mode is reached when fc,. = 0. We also note, that in the excitation of both It- and t-wave it is 
the X component of the electric field that is growing exponentially. 

Estimating and Q using ^-function ( max[|^] « 0^25 /ua and max[(^^^|^ " T^^)] ~ v^)' 
we wind the maximum growth rates of the t- and It-modes in the limit 5 3> l/7^g^: 

H..res5 ( 7^ 9/(7) A 

.i + "d7Vc2 97 



where Wp res is the plasma density of the resonant particles. 



We estimate the growth rates (B7) for the distribution function of the resonant particles having a 
Gaussian form: 

where pb is the momentum of the bulk motion of the beam and Ap is the dispersion of the momentum. 
Assuming in (p7| ) that ua^b /c 1 we find the growth rates 

r-r'-yi%^4, (B9) 

where 7^ = pi,/{mc), A7 = Ap/lmc). 



C. Conditions on Cherenkov-drift instability 

Since Cherenkov-drift resonance requires a very high parallel momentum, the resonant interaction 
will occur on the high phase velocity waves. This implies that, like cyclotron-Cherenkov resonance, 
Cherenkov-drift resonance is always important for the extraordinary mode, while for ordinary and Alfven 
modes it is important only for small angles of propagation. We can then use the small approximation to 
dispersion relations for the small angles of propagation ( |Al[ ) . Introducing cylindrical coordinates x ,r 
with r along the radius of curvature of the field line, x perpendicular to the plane of the curved field line 
and (f) the azimuthal coordinate (Fig. the resonance condition (||) may then be written as 

^ S+U' + Uo-Ud/cf = (CI) 



272 " ' 2'" '2 

Ires 

where we used Vres = c(l — ^), ^ = and ^ — 



We note that for the extraordinary mode S is independent of frequency and equation (CI) becomes 
independent of frequency as well. This means that for the given angle of propagation Cherenkov-drift 
resonance on the extraordinary mode occurs at all frequencies simultaneously. It is different from the 
cyclotron-Cherenkov resonance occurs at a particular frequency. For both ordinary and Alfven modes 6 
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does independent on the frequency, so that for a given angle of propagation the Cherenkov-drift resonance 
occurs at a fixed frequency. 

Let us now discuss the condition for the development of the Cherenkov-drift instability. First we drop 
the — term from equation ( |Cl| ). This term is much smaller than 6 for the radii satisfying 



R 



> 



* 3 



1/3 



which is satisfied everywhere inside the pulsar magnetosphere. 



We find then that Eq. (CI) can be satisfied for 

^ < V2d, \9 - Ud/c\ < V2S 
Fig. 1^ describes the emission geometry of the Cherenkov-drift instability. 



From (C3) we see that drift of the resonant particles becomes important if 



Using the expression for S (M), Eq. (C4) gives 



\RnsJ 



7p lb 



(C2) 



(C3) 



(C4) 



(C5) 



At a given radius this condition may be considered as an upper limit on the curvature of the field lines 

(C6) 



2 3 2 \ 

Rb < I ^S^Z I « 2 X lO^cmfori? = lOOOi^^vs 



Alternatively, condition ( |C5| ) may be regarded as a lower limit on the radius from the star. In the 
dipole geometry the radius of curvature for the open field lines may be estimated as Rb > \/ ^ The we 



find from Eq. (|C5D 



\RnsJ \ C7372 J 



(C7) 



In what follows we assume that the Cherenkov-drift resonance occurs in the outer parts of the pulsar 
magnetosphere where the typical value of the drift velocity is Ud « 0.01c. 

The growth rate for the Cherenkov-drift resonance instability on the primary beam is calculated in 
appendix ^ Eq. B9. Expressing the growth rate BE in terms of the parameters of the magnetospheric 
plasma we find 

^ AT„c^ ^^^^ 



The growth rate (C8) should satisfy several conditions. The first is the criteria for the fast growth: 
T/il > 1. This is the requirement that the growth rate is fast enough, so that the instability can have time 
to develop before the plasma is carried out of the magnetosphere. From ( |C^ ) we find that the condition 
r/rj > 1 is satisfied for the chosen parameters and the typical frequency of emission uj — 5 x 10^ rad s~^. 
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The next condition that a growth rate should satisfy is that the growth length be much less that the 
length of the coherent wave-particle interaction (A8). Estimating the range of resonant angles sa \/^ 



and using the growth rate (|C8|), this condition gives 



Rb > 



lO^^cm at R ' 



10 cm 



(C9) 



Equation (C9) is the lower limit on the radius of curvature of the field lines. It will restrict the emission 
region to the field lines closer to the central field line, where the radius curvature is large. 



There is a third condition that the growth rate (C8) should satisfy - the condition of the kinetic 
approximation (All). For the Cherenkov-drift resonance condition (All) gives 



kxUd 

7 



>r 



(CIO) 



Estimating k^c ~ uJUd/c we find that this condition can be easily satisfied in the pulsar magnetosphere. 

In this Appendix we showed that the Cherenkov-drift instability can develop in the pulsar 
magnetosphere. The region of the development of the instability has a radius of curvature limited both 



from below (by the condition of a large radius of curvature C6) and from above (by the condition of a short 
growth length |C9|). 



D. Instabilities in millisecond pulsars 

Both cyclotron- Cherenkov and Cherenkov-drift instabilities may develop in millisecond pulsar. It 
is harder to satisfy the conditions for the development of the cyclotron-Chcrenkov instability in the 
millisecond pulsars than in the normal pulsars. Since there is no clear distinction between core and 
cone-type emission for the millisecond pulsars it is possible that only Cherenkov-drift instability develops in 
their magnetospheres. Alternatively, different pulsars may have a substantially different plasma parameters, 
so that both instabilities may develop in some of them. 

Here we will discuss briefly the conditions for the development of thes instabilities in a "standard" 
millisecond pulsar with the period P — 5 x 10^^ s and the surface magnetic field B = 10® G. As a first 
order approximation we will assume that the other plasma parameters, i.e. the initial primary beam 
Lorentz factor 7^"'' = 6 x 10^, the primary beam Lorentz factor at the light cylinder 7^"'' = 6 x 10^, 
its scatter A7 = lO'^, average streaming energy of the secondary plasma 7^ = 10, its scatter in energy 
Tp = 10 and the multiplicity factor A = 3 x 10^ are the same. The light cylinder is now at a radius 
Ric = 2.4 X lO^cm = 2ARns- These are very approximate assumptions. The plasma parameters in the 
millisecond pulsars are likely to be different from the normal pulsars. At this moment we do not understand 
the physics of the pair production well enough to make a quantitative distinction between normal and 
millisecond pulsars. 

The growth rate of the cyclotron-Cherenkov instability becomes equal to the rotational frequency of 



the pulsar at (see ( A6)) 



which is inside the light cylinder. 
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The approximation of the kinetic growth rate cyclotron-Cherenkov instability (A13) requires that 



\Rns J 



< 3 X lO'' 



(D2) 



which is satisfied. The condition of a short growth length of the cyclotron-Cherenkov instability (A8) 
requires that 



Rb > lO^cm 



(D3) 



at the emission site of i? « 20 Rns- This condition is hard to satisfy: cyclotron-Cherenkov instability does 
not develop in the millisecond pulsars. 



For the Cherenkov-drift instability, the condition of a large drift (C4) is now satisfied for the all open 
magnetic field lines as well as the conditions of a fast growth (condition (i) on page |lj) ans the condition 
of the kinetic approximation (condition (ii) on page Cherenkov-drift instability does develops in the 
magnetospheres of the millisecond pulsars. 



E. Effects of curvature radiation reaction on the beam 



Consider a beam of electrons propagating in s dipole curved magnetic field. The radius of curvature on 
a field line near the the magnetic moment is 



4 y/RNsR 



(El) 



Here a* ^ ^JRns/R is the angle at which a given field lines intersects the neutron star surface. 

An evolution of the distribution function under the infiuence of the radiation reaction form is described 
by the equation: 

dfiz,pz,t) d (dpz \ 



Equation (E2) can be solved by integrating along characteristics 

26274^3 



dt 



ic'Rl 



(E3) 



With the radius of curvature given by (El) Equation ( |E3| ) has a solution 



1 / \/p1 + {mcY 



In 



1 + v^f+T^ 



Po 



3 r^a 



*2 



2 Rns 



Ht/to) 



(E4) 



where = is the classical radius of an electron, Iq ~ Rns/c and po is the initial momentum. This 
equation can be solved for pz 3> mc 



-1/3 



(E5) 
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Solution of the continuity equation (E2) is then 

/(^,7,t) = /o(^7o(7,i), ^ = *") 

where foi^Oit — to) is the initial distribution function at i? = Rns- Evolution of the distribution function 
is shown in Fig. 

If, originally, the beam had a scatter in energy A70 then at a later moment it will have a scatter in 
energy given by 

A7 = A70 [I + ^ojj^ In(V^o) j (E7) 

Which implies that the scatter in energies of the beam may decrease much faster than its average 
energy. 

To estimate the decrease in the average energy and in the energy scatter onsider an evolution of the 
beam on the last open field line. Then a*"^ ~ Rns^/c and at the light cylinder (t = 2Tr/fl) we have 



7 f ^ ^ 39 rgfi J ^ 2TTC 



-1/3 



Which may be expressed as 



A70 V 2 c 'Rns^ 

A7 A70 / 7 



-4/3 



A7 / ^9refl , , 27rc , , , , 

= 1+ 7o7T— li^(l^) (E8) 



(E9) 



7 70 V7o, 

If the original beam was mildly relativistic ( A70 « 0.1 70 then in order to reach a scatter A7/7 w 10^^ 
the beam have to lose about 90% of its original energy: 7/70 ~ 0.1. Estimating ( p^ ) with we find that the 
scatter in the energy of the beam on the last open field line at the light cylinder for the normal pulsars may 
be as small as A7ic = 100. This vindicates our assumption of the chosen beam energy spread. 

In the normal pulsars the cooling of the high energy beam is important for the 70 > 2 x 10^. In 
millisecond pulsars it becomes important for 70 > 2 x 10^. The cooling of the resonant particles increases 
the growth rates of both cyclotron- Cherenkov and Cherenkov-drift instabilities. On the other hand, the 
kinetic approximation condition for the Cherenkov-drift instability may not be satisfied for a very small 
scatter in energy of the resonant particles. Unlike the cooling, the change in the average energy of the beam 
is normally not very important and is neglected. 



F. Instability in an ion beam 

If the relative orientation of rotational and magnetic axis is such that the electric field pulls up positive 
charges, then the primary beam will consist of ions. The critical Lorentz factor of the charge particles 
needed for the pair production is independent of the mass of the particle. Ions will be accelerated to the 
same Lorents factors as electrons 7^ « 10^^^. It is natural to expect that in approximately 50% of neutron 
stars the primary beam will consist of ions. In what follows, the ratio of ion and electron masses will be 
denoted as p = rrii/me. 
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First we consider a cyclotron- Cherenkov instability on the ion beam. If the primary beam consists 
of ions, then the cyclotron-Cherenkov resonance at the given frequency occurs closer by a factor p^/^ w 3 



(see A4). This will increase the density of the resonant particles by ~ 30, but the growth rate for the 



cyclotron-Cherenkov instability (A5) is proportional to the squared plasma frequency of the resonant 
particles, which is inversely proportional to the mass of the particles. Taking into account that the condition 
of the kinetic instability is independent of the mass of the particle, we conclude that the cyclotron-Cherenkov 
instability on the ion beam does not develop due to the very small growth rate. 



The growth rate of the Cherenkov-drift instability is proportional to p"^ (Eq. B£), so the growth rate 
on the ion beam is many orders of magnitude smaller than the growth rate on the electron beam. Thus, we 
conclude that both cyclotron-Cherenkov and Cherenkov-drift instabilities do not develop on the ion beam. 
The pulsars with such a relative orientation of magnetic and rotational axis that the induced electric field 
accelerates ions can produce radio emission only at the cyclotron-Cherenkov instability on the tail particles. 
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Fig. 2.— 



Example of emission geometry that produces core and cone profiles. 
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Fig. 3. — The dotted lines show constant spin down ages and magnetic fields. There is also a dashed line 
for the equilibrium period spin-up. The heavy dashed curves and shading illustrate several death lines from 
the literature (see Hansen 1997 for references). The circles indicate binary pulsars and those points with 
horizontal lines through them are in globular clusters with apparently negative E (they are contaminated 
by cluster accelerations). 
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Fig. 4. — Chcrenkov-drift emission. Drift velocity is perpedicular to the plane of the curved field line 
(B — Rc) plane, Rc is a local radius of curvature). The emitted electromagnetic waves are polarized along 
Ud- The emission is generated in the cone centered at the angle 9^"^ = Ud/c and with the opening angle 
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Fig. 5. — Dispersion curves for the waves in a cold electron-positron plasma in the plasma frame in the limit 
LOp <^ LOB- There are three modes represented by the dashed (ordinary mode), solid (extraordinary mode) 
and long dashed (Alfven mode). The dotted line represents the vacuum dispersion relation. For the exact 
parallel propagation, the dispersion curves for the ordinary mode and Alfven mode intersect. The insert in 
the upper left corner shows the region near the cross-over point wo- 
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Fig. 6. — Coordinate system. 
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Fig. 7.— 



Emission geometry of the Cherenkov-drift resonance 
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Fig. 8. — Top view of the field lines in the equatorial plane of a rotating point dipole in vacuum ( Yadigaroglu 



1997). Circle indicates light cylinder. The locations of cyclotron-Chercnkov and drift instabilities are shown 



(similar regions will be on the other side of the pulsar). The cyclotron- Cherenkov instability develops in the 
region of almost straight field lines. The location of the Cherenkov-drift emission depends sensitively on the 
curvature of magnetic field line. Two possible locations of the Cherenkov-drift emission are shown: ringlike 
near the magnetic axis and in the region of swept field lines (shaded ellipse). When the effects of plasma 
loading are taken into account, the field lines will become more curved. The current flowing along the open 
field lines will also produce a torsion in the field lines, so that at the light cylinder the structure of the field 
will be changed considerably 
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Fig. 9. — Switch of the sense of the circular polarization due to the cyclotron- Cherenkov resonance on the 
tail of the plasma distribution. The drift velocity of the electron and positron is in the opposite direction. 
As the line of sight crosses the emission region, the observer first will see the waves emitted by the particles 
of one sign of charge and then the other. The waves emitted by electron and positron have different circular 
polarization. 
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Fig. 10. — Evolution of the primary beam distribution function due the curvature radiation reaction. For 
the illustrative purposes we chose the initial beam 7-factor equal to 6 x lO'', and the initial scatter in Lorentz 
factors A7 = lO'^. 



